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ONE  NEAR-TERM  RESPONSE  to  current  constraints  of 
fuel  cost,  emissions  standards,  and  acceptable 
performance  is  near-lean-limit  operation  of  con- 
ventional engines.  However,  the  reduced  pressure 
and  temperature  achieved  in  homogeneous-charge 
four-stroke  spark-ignition  engines  under  such 
fuel-lean  operation  may  lead  to  enhanced  emission 
of  partially  oxidated,  or  pyrolyzed  but  unburned, 
hydrocarbons  known  to  be  toxic  and/or  carcino- 
genic [1].*  Proper  design  of  the  number,  loca- 
tion, timing,  and  energy  of  spark-ignition 
sources  should  preclude  bulk-gas  flame-out  as  a 
source  of  such  emissions.  However,  the  quench 
layers  formed  on  cold  cylinder  surfaces  (head 
wall,  side  wall,  and  piston  crown)  [2]  are 
thickened  under  off-stoichiometric  operation  [3] ; 
unless  the  cold  unburned- hydrocarbon  content  of 
these  quench  layers  can  be  oxidated  by  mixing 
with  hot  oxygen-rich  bulk-gas  products  before 
being  exhausted,  the  lean-operation  strategy  may 
be  flawed.  What  seems  required  are  experimental 
and  theoretical  investigations  of  the  formation 


The  unsteady  one-dimensional  interaction 
of  a planar  flame,  propagating  through  a fuel- 
lean  premixture,  with  a parallel  impervious 
noncata lytic  wall  is  modeled  as  a Stefan-type 
problem  Incorporating  one-step  chemistry. 
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and  fate  of  unburned  hydrocarbons  in  reciprocal ng- 
piston-type  internal  combustion  engines  [4],  to 
furnish  as  much  engineering  information  as  pos- 
sible for  the  new  designer  who  faces  decisions 
concerning  the  relative  merits  of  alternative 
engine  modifications.  This  paper  is  one  in  a 
series  in  which  the  authors  seek  to  contribute 
to  this  need. 

The  scope  of  the  pertinent  investigation 
should  be  broadened  in  view  of  recent  interest 
in  low-heat-transfer  zircomium-oxide-sprayed 
engine  components,  and  In  adiabatic  components 
of  glass  ceramics  (such  as  silicon  nitride  and 
silicon  carbide),  particularly  for  larger 
engines  [5,6].  If  durable  materials  can  be 
developed,  not  only  are  problems  associated  with 
the  size  and  weight  and  reliability  and  mainte- 
nance of  the  cooling  system  reduced  or  even 
eliminated,  but  also  currently  wasted  heat  can 


*Numburs  in  brackets  designate  references 
given  at  the  end  of  the  paper. 


acceleration  of  flame  speed  near  an  adiabatic 
wall  (pertinent  to  proposed  engines),  are 
obtained.  The  time  history  of  unburned  hydro- 
carbon content  varies  with  dimensionless 
Arrhenius  activation  energy,  the  order  of  the 
reaction  chemistry,  and  the  Lewis-Semenov 
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be  utilized  via  engine  compounding  (if  the 
exhaust  port  of  the  primary  engine  can  survive 
hotter  operating  conditions). 

Thus,  interest  centers  on  interaction  of 
flames  propagating  through  confined  fuel -lean 
premixtures  with  (l)  isothermal  (cold)  walls, 
and  (2)  adiabatic  walls.  Burning  is  anticipated 
to  be  retarded  near  cold  walls  and  accelerated 
near  adiabatic  walls,  but  quantitative  insight 
relating  results  to  controllable  parameters  is 
extremely  limited  in  the  existing  literature. 

It  seems  premature  to  undertake  the  turbu- 
lent case  when  results  for  the  laminar  case  are 
not  yet  available.  Further,  formulation  for 
turbulent  reacting  flow  in  an  automotive  context 
is  today  highly  conjectural  [7],  and  turbulent 
fluctuations  may  be  significantly  reduced  in  the 
near-wall  phenomena  of  interest.  In  any  case, 
interest  here  is  centered  on  laminar  flow, 
though  unsteady  and  nonisobaric  conditions  are 
ultimately  to  be  included. 

Though  intermediate  orientations  can  and  do 
arise  in  practice,  it  is  convenient  to  confine 
attention  to  two  arrangements  [8,9]:  one  in 
which  the  bulk-gas  isotherms  are  parallel  to  the 
containing  wall  (henceforth  termed  the  end-wall 
geometry)  [10-14],  and  one  in  which  the  bulk-gas 
isotherms  are  perpendicular  to  the  containing 
wall  (henceforth  termed  the  side-wall  geometry) 
[15,16].  Elsewhere  the  authors  intend  to  present 
both  theoretical  and  experimental  results  for  the 
side-wall  geometry.  It  appears  appropriate  to 
confine  the  present  discussions  to  theoretical 
results  for  the  end-wall  geometry  only. 

In  general,  attention  is  confined  to  un- 
steady one-dimensional  flow  in  which  a well- 
developed  flame  initially  is  propagating  through 
a fuel -lean  premixture;  both  the  nonconfined 
(isobaric)  and  confined  (nonisobaric)  cases  are 
oi  interest.  How  the  properties  of  the  flame 
become  modified  as  a wall  is  approached  is  the 
central  issue.  Nevertheless,  as  a prerequisite, 
some  properties  of  flame  propagation  away  from 
walls  must  be  established  to  serve  as  the  start- 
ing (and  reference)  conditions.  The  dynamics 
is  not  central,  once  it  is  noted  that  the  pres- 
sure is  spatially  invariant,  though  not  in 
general  temporally  invariant,  for  the  highly 
subsonic  phenomena  of  interest.  What  is  central, 
in  addition  to  overall  conservation  of  mass,  is 
conservation  of  species  and  of  heat;  this 
conservation  consists  of  a balance  of  convection, 
diffusion,  and  chemical  reaction. 

A conventional  set  of  self-consistent 
approximations,  often  alluded  to  (in  toto)  as  the 
Shvab-Zeldovich  formulation  [17],  is  adopted; 
this  formulation  has  been  of  great  engineering 
value  in  the  analysis  of  a wide  variety  of  aero- 
thermochemical  phenomena.  A mixture  of  ideal 
gases  of  comparable  molecular  weight  is  taken  to 
be  present,  such  that  a universal  binary  diffu- 
sion coefficient  suffices,  and  one  (constant, 
order-unity)  Lewis-Semenov  number  may  be  intro- 


duced to  characterize  the  ratio  of  thermal 
diffusivity  to  mass-transfer  coefficients. 

The  thermal  diffusivity  is  taken  to  be  inversely 
proportional  to  the  square  of  the  density.  A 
universal  constant  heat  capacity  suffices,  and 
thermodiffusion,  barodiffusion , radiative  trans- 
fer, and  mechanical  dissipation  of  energy  are 
negligible.  Most  notably,  a direct  one-step 
irreversible  exothermic  mechanism  describes  the 
chemical -reaction  mechanism,  with  an  Arrhenius 
form  appropriate  for  the  specific  rate  constant 
occurring  in  the  law  of  mass  action  [3],  The 
fact  is  that  the  detailed  rates  and  multistep 
mechanisms  of  the  actual  higher-hydrocarbon/ait 
chemical  reaction,  particularly  with  hetero- 
geneous complications  posed  by  the  presence  of  a 
wall  [18],  are  not  known;  once  this  chemical- 
kinetic  compromise  is  accepted  as  currently 
necessary,  the  other  simplifications  follow  from 
the  considerations  of  a consistent  level  of 
approximation. 

Accordingly,  aside  from  the  role  of  chemical 
reaction,  which  introduces  the  complexity  of 
transcendental  nonlinearity  into  the  pertinent 
initial/boundary-value  problems,  treatment  of 
parabolic  linear  constant-coefficient  equations 
only,  familiar  from  the  field  of  heat  conduction, 
is  required.  However,  for  Arrhenius  activation 
temperature  large  relative  to  burned-gas  tempera- 
ture, the  role  of  chemical  reaction  is  limited  to 
a spatial  domain  small  relative  to  the  domain  in 
which  a chemically  frozen  balance  of  convection 
and  diffusion  is  an  adequate  description  [19-21]. 
Here,  this  small  spatial  domain  is  reduced  to  a 
mathematical  interface  [22],  which  is  a Dirac- 
del ta-function-type  sink  for  reactants,  and 
source  for  products  and  chemical  heat  release. 

The  interface  is  then  a premixed  thin  flame, 
analogous  to  the  Burke-Schumann  diffusion  flame 
[23]  long  employed  in  aerothermochemical  flows 
involving  unpremixed  reactants;  however,  the 
flame  without  structure  in  the  premixed  case,  in 
general,  retains  chemical-kinetic  parameters 
(such  as  the  Arrhenius  activation  energy),  while 
the  thin  diffusion  flame  retains  no  chemical- 
kinetic  parameters.  Only  in  the  special  case  of 
an  adiabatic  flow  at  Lewis-Semenov  number  unity 
is  the  premixed  flame  without  structure  indepen- 
dent of  chemical -kinetic  parameters  [13]. 
Nevertheless,  adoption  of  a thin  flame  (at  which 
temperature  and  mass  fractions  are  continuous, 
but  their  derivatives  are  discontinuous)  is  a 
significant  simplification,  even  for  the  premixed 
reacting  flow,  because  an  overall  solution  is 
obtained  simply  by  appropriately  joining  local 
solutions  to  heat-conduction-type  problems  at 
interfaces  (flames).  In  that  the  position,  as 
well  as  the  temperature,  of  the  flames  must  be 
found  in  the  course  of  solution,  thin-flame 
modeling  reduces  intractable  aerothermochemical 
problems  to  tractable  Stefan  problems  [24,25]. 

Not  only  are  the  conventional  analytic  proce- 
dures for  treating  Stefan  problems  available, 
but  the  authors  have  introduced  novel  techniques 
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for  the  aerothermochemical  applications  of  rele- 
vance here.  Such  modeling  may  prove  useful 
in  a wider  range  of  premixed  reacting  flows  than 
the  particular  flame/wall  interactions  of  in- 
terest here. 

FORMULATION 

Consider  the  unsteady,  one-dimensional, 
low-Mach-number  flow  of  a reacting  premixture, 
where  the  chemical  reaction  is  a direct,  one- 
step,  irreversible  one  between  oxidant  0 and 
fuel  F that  generates  gaseous  product  P,  namely: 

v0°**FF-vpP.  (!) 

' where  is  the  stoichiometric  coefficient  of 

species  i,  i = 0,  F,  P.  In  the  Shvab-Zeldovich 
approximation,  the  appropriate  dimensional 
equations  of  conservation  of  mass,  conservation 
of  species,  conservation  of  energy,  and  of 
state,  for  such  a flow,  are: 


(2) 


■(»-&)•*(--&)■ 


°*(f£ 4 u*  &)-  & !£)*  (q*/c*pw*  * ; <4) 

p*  - c*  R*  T*.  (5) 

Here,  t*  and  x*  are  the  time  and  Cartesian 
spatial  coordinates,  respectively;  u*  is  the 
gas  speed  and  y*  is  the  streamfunction;  p*  and 
T*  are  the  density  and  temperature;  R*  is  the 
gas  constant  for  the  mixture,  taken  to  be  com- 
posed of  species  of  comparable  molecular  weight; 
p*  is  the  pressure,  taken  to  be  a function  of 
time  only,  from  consideration  of  the  equation  of 
conservation  of  momentum;  Y.  is  the  stoichio- 

metrically-adjusted  mass  fraction  of  species  i; 
D*  and  X*  are  the  mass-transfer  and  thermal - 
conductivity  coefficients,  respectively,  and  c* 
is  the  (constant  universal  ) heat  capacity  at  p 
constant  pressure;  and  Q*  is  the  specific  heat 
of  combustion;  while  w*  is  the  reaction  rate. 

In  what  follows,  it  is  taken  that  the  Lewis- 
Semenov  number  is  a constant  of  order  unity, 
i.e. , 


Le  ■ * const-  * 00).  (6) 


c*2  0 « const.  ( 7 ) 

In  this  presentation,  the  law  of  mass  action, 
under  an  Arrhenius  specific  rate  constant, 
modified  to  overcome  the  so-called  cold-boundary 
difficulty,  is 

«*  * B-  T*°  c*V°  "F  YfVf  Y0V°  exp  j-0*/(T*  - TJ)j  . (6) 

where  B1*  is  the  (constant)  frequency  factor; 

0*  is  the  (constant)  Arrhenius  activation 
temperature;  T*  is  the  cold-boundary  (i.e.,  cold- 

premixture)  temperature;  and  a characterizes  the 
pre-exponential  thermal  dependence  of  the 
reaction  rate.  While,  in  general,  there  need  be 
no  association  of  this  dependence  and  order  of 
reaction,  for  convenience  (and  within  the  accu- 
racy of  previously  adopted  approximations),  here 
it  is  taken  that 


« * (v0  * VF>  * '• 


(9) 


It  is  also  convenient  to  introduce  the 
von  Mises  transformation  (x*,t*)  -*■  (>l^,t*), 
where  the  streamfunction  t|>*  is  defined  in  Eq.  2. 
Under  this  transformation,  the  species  and 
energy  conservation  equations.  3 and  4,  become 


st*  - (o*2°*)  ~7  * -B'*(p*/r*)<v°*  F>"'  r/r  Y0V°  exp 


!-w  • 4(iii 


£ - <‘-’n  '•  ~i  ■ "«;>  g 


c»)  B'»(pVR*) 


(V0*VF>’1 


-0*/{T*  - T*)J 


while  the  (complementary)  mass  conservation  and 
mapping  equations  are  given  by 


9y*  , » /_L\  . >»l . ' . 

W at*  \o*)  * a**  c* 


(12a, b) 


To  develop  a nondimens ional  formulation,  let 


Further  it  is  compatible  with  the  accuracy 
characteristic  of  the  Shvab-Zeldovich  approxi- 
mation to  adopt 
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(16) 


Here,  the  subscript  u denotes  unburned  (or  cold- 
premixture)  conditions;  while,  the  subscript  b 
denotes  burned  conditions  (in  the  absence  of 
heat  losses,  i.e.,  under  adiabatic  conditions). 
More  specifically,  u*  is  the  so-called  adiabatic 

flame  speed.  Also,  * = (YFu/Y0u)  is  the  (so- 

called)  equivalence  ratio,  which  is  less  than 
unity. for  cases  of  interest  here.  In  this 
development,  it  is  taken  that  Ypb  = 0,  such  that 

Y0b  = (Y0u  " YfJ  = t(l  - 0)/<t>3Ypu , and  that 
Tb  = + (Q*/c£)YFu].  Further»  = R*  T*. 

Thus,  the  nondimens ional  species  and  energy  equa- 
tions can  be  written  as 


9 ' m * » 

» * p 


where 


, ill  0 ♦ m)  dg 

a <5  dt 


a;  - (o'/cjjr, 
K . —S.  " 1 


(18) 


6 * - TJl  * (20) 

'(^)(s)vv'’or>:rv"-..-  («) 


steady  adiabatic  laminar  flame  propagation  for 
large  6,  it  is  known  (19-21]  that  (1)  the  flow 
field  may  be  divided  into  a chemically  frozen 
preheat  domain,  of  scale  (D*/u*),  in  which  con- 
vective diffusion  occurs,  and  a relatively 
thinner  reactive  domain,  of  scale  6_1(D*/u*),  in 

which  chemical  reaction  and  diffusion  balance; 

(2)  that  (except  for  very-near-stoichiometric 
conditions,  i.e.,  <j>  ■*  1 ) , 

'♦tV*'  (23) 


in  the  reaction  zone,  the  only  region  in  which 
chemical  processes  play  a significant  role  (so 

Y0  ~ Y0  for  a11  'M  from  Ecl-  14);  and  (3)  that 
u 

(if  r(z)  denotes  the  (complete)  gamna  function  of 


(24) 


for  the  fuel-lean  case.  Accordingly,  for  more 
complicated  phenomena  involving  combustion  in 
fuel-lean  premixtures  characterized  by  large  6, 
it  appears  convenient  and  appropriate  to  reduce 
the  narrow  reaction  zone  to  a mathematical  inter- 
face at  iji  = H* ( t ) , with  H* ( t ) , the  interface 
location,  a quantity  to  be  found  in  the  course  of 
solution.  For  continuity  of  the  dependent  varia- 
bles themselves,  it  is  required  that 


y - o,  t - R ■ R(t)  as » ->  *(t),  for  t ; o , (25) 


with  H(t),  the  temperature  at  the  flame,  another 
quantity  to  be  found  in  the  course  of  solution. 
The  flame  is  anticipated  to  be  a Dirac-delta- 
function-type  source  [26]  for  chemical  heat  and 
sink  for  reactants,  such  that  first  (and  higher) 
spatial  derivatives  of  the  dependent  variables 
(aside  from  the  pressure)  are  discontinuous. 
Hence,  in  the  vicinity  of  t(t),  the  equations  for 
conservation  of  fuel  and  energy  are  well 
approximated  by 


The  nondimens ional  state,  continuity,  and  mapping 
equations  can  be  written  as 


^ P 


♦ w)  - 1 


*}• 


(26) 


U 4 • -A  P 
* 


r>  - ' 


(27) 


Typically,  the  nondimensional  activation 
temperature,  B,  is  much  greater  than  unity  [3]; 
this  observation  is  central  to  the  analysis  to 
follow.  The  quantity  A plays  the  role  of  an 
eigenvalue  in  steady  flame  propagation;  again, 
it  is  steady  adiabatic  flame  propagation  into  a 
quiescent  premixture  characterized  by  p*,  T*, 

YFu*  and  that  yields  uj,  taken  as  known  for 
species  0,  F reacting  to  form  product  P.  For 


Hence,  in  the  vicinity  of  r(t), 

[1 . lt  T)  . 0 — > iJilLLlilal . aiiOilaJL  »natkt2  (28) 
a,.*  | H h | a* 

if  the  convention  is  henceforth  adopted  that 
Y_[^>  *y(t)  ,t  1 = OT  It  suffices  for  the  calc  u 1 a - 
tion  of  the  dependence  of  the  jump  in  Y^  on  B and 

H(t)  (using  Eq.  26)  to  take 
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T <•»*)  • H(«)  • »(*.«)  ♦ ...  . (29) 

The  discrepancy  is  of  the  order  of  1/u  and  would 
be  eliminated  formally  in  a derivation  using 
singular  perturbation  theory.  Accordingly,  for 

e <<  1 1 

f(t)H  , . , . » , »(t)H 

*(«).«  * fit)..  * * 

or.  from  Eq.  29. 

. n „(vo  * • »/  wj.,  IMLilIAilj  (31 } 

where  the  upper  limit  has  been  extended  to  infi- 
nity by  use  of  Laplace-integral-type  (i.e., 
steepest-descent)  concepts  [261.  Expansion  of 
the  argument  of  the  exponential  gives 

M’ — ’'‘v' "•+¥!/■ 

: Jlv *V  ,J.  I Hj . (33) 

by  Eq.  24. 

These  results  may  be  anticipated  to  hold 
for  vigorous  burning  only,  such  that  a narrow 
reaction  zone  may  be  adequately  approximated  by 
equivalent  interfacial  conditions.  Nevertheless, 
the  results  should  be  of  qualitative  use  for 
less  vigorous  burning,  in  which  the  reaction 
zone  is  diffuse.  The  interfacial  results  pre- 
sented here  by  heuristic  arugment  may  be  obtained 
more  formally  by  application  of  multiple-scaling 
techniques  of  modern  asymptotic  analysis; 
actually,  Zeldovich  and  Spalding  intuitively 
used  equivalent  methodology  decades  ago  [17]. 

The  resulting  Stefan-type  problem  [24,25] 
may  be  summarized  as  follows  for  the  domain 
t > 0,  u»L  < ip  < n,R.  For  t(t)  < 4*  < 4>R,  i.e., 

downstream  of  the  flame,  the  species  and  energy 
equations  are  (Fig.  1) 


T-l  I ♦ t M 

1 T p ( 


while  for  4^  < 41  < Y(t),  I.e.,  upstream  of  the 
flane,  the  species  and  energy  equations  are 


S-9-*- 


>H  1.  y N . Ill  \ * *H  dp 

n ’ S?  ' at 


u ■ S • 0 


T . 0 or  g ■ 0 


premixture 
1 > Y'  0,  H >T  >0 


flame:  Y * 0,  T = H(t) 


4-> 

%-t 
] 

-*-> 

— 1 ♦ 

burned 

gas 

Y = 0 

xf(t) 

x x 

¥ * “* 

¥(t) 

» 

0“1 

• 0 

0 

Fig.  1 - Schematic  for  a thin-flame,  Stefan-type 
model  of  the  unsteady  one-dimensional  isobaric 
flow  in  which  a large-activation-energy  laminar 
flame  propagates  toward  an  impervious  noncata- 
lytic  wall,  that  is  either  isothermal  (cold)  or 
adiabatic.  The  purely  diffusive-convective 
unburned-premixture  region  is  characterized  by 
finite  fuel  mass  fraction  Y and  by  temperature  T 
below^the  temporally  varying  thin- flane  tempera- 
ture H(t);  the  purely  diffusive-convective  burned- 
gas  region  possesses  no  fuel.  In  the  dimension- 
less Cartesian  coordinate  x,  the  flame  site  is 
denoted  Xf ( t ) ; in  the  von  Mises  (mass-weighted) 
coordinate,  the  corresponding  variables  are  4s 
Y;  in  the  Landau  (flame-stabilized)  coordinate, 
the  corresponding  variables  are  0,  1.  The  gas 
speed  at  the  flame  is  denoted  uf;  the  flame  speed 
(in  laboratory  coordinates)  is  denoted 


The  following  notation  has  been  Introduced  for 
convenience: 


Ht*t)  for  » * t(t) 


for  , * t(t). 


At  the  interface,  i.e.,  for  4’  ♦ Y(t)  (t  ' 0), 
the  boundary  conditions  for  the  dependent  varia- 
bles are 

» - 0.  ».m  . R(t)  ( 35 

while  the  boundary  conditions  for  the  first 
derivatives  of  the  dependent  variables  are 


[(s  * l*  5?)  • u S],  * 


(<v0  ♦ vF)  - i)/:  v,  ♦ 


’ 4 * °?) 
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For  t * 0,  the  initial  conditions  are  [with 
(to  * such  that  the  f^a,ne  is  initially 

far  from  the  unburned-premixture  boundary] 


¥ * »q.  Qiven,  R • Ng,  given,  p * 1. 


(42) 


where  bp,  Cq,  dQ  are  given  consts.  close  to 

unity  in  value,  so  the  initial  profiles  are  very 
similar  to  those  occurring  in  laminar  isobaric 
flame  propagation.  The  right  end-point  boundary 
condition  is 

h ♦ i « v - *,<t : C);  (45) 

while  the  left  end-point  boundary  conditions  are 


- o »s  * - ,t(t  ; 0). 

(46) 

U - o « * - <L(t  ; o). 

(47) 

For  some  cases  to  be  considered  below,  the  boun- 
dary and/or  initial  conditions  given  here  are 
modified,  but  explicitness  is  preferred  to  gene- 
rality for  purposes  of  presentation.  However, 
any  initial  conditions  adopted  are  anticipated 
to  be  compatible  with  the  boundary  conditions  to 
a high  degree  of  accuracy.  For  brevity,  the  com- 
plementary relations  (Eq.  22)  are  not  repeated 
here,  though  it  is  noted  that  boundary  conditions 
on  the  velocity  at  impervious  walls  are  enforced, 
of  course,  in  the  obtaining  of  solutions. 

One  special  case  is  deemed  worthy  of  dis- 
cussion at  this  point.  For  a semi-infinite 
domain  (e.g.,  -*■»),  the  flow  is  well  approxi- 

mated as  isobaric  in  many  circumstances,  such 
that  p(t)  = 1.  If,  in  addition,  the  flow  is 
adiaKcic  with  Lewis-Semenov  number  Le  = I,  then, 
under  the  provision  that  the  following  statements 
are  compatible  with  the  initial  and  boundary 
conditions  of  interest. 


h • H ■ 1 for  » ; T(t);  Y ♦ H • R • I for  . ; H . 


In  this  special  circumstance, 


-l  — > - 1 « » ■>  »\ 


(48) 


(49) 


independently  of  the  chemical-kinetic  parameter 
6.  Because  of  this  significant  simplification, 
this  special  adiabatic  case  has  been  considered 
separately  for  the  particular  circumstance  of  a 
flame  propagating  isobarically  toward  an  adia- 
batic, noncatalytic  wall  [13]. 

In  the  sense  that  the  reaction  order  and 
Arrhenius  activation  temperature  for  a one-step 
pseudomechanism  for  the  chemistry  must  be  empi- 
rically assigned,  the  exponent  on  the  pressure 
factor,  the  exponent  on  the  pre-exponeptial 
temperature  factor,  and  the  nondimensional  acti- 
vation temperature  in  the  exponential  factor, 
all  in  the  expression  for  the  mass-fraction 
gradient  at  the  flame,  are  available  for  assign- 
ment [3,27]. 

Solution  to  the  Stefan-problem  formulation 
may  be  sought  by  direct  numerical  integration; 
this  can  be  a challenging  task  to  accomplish. 
Thus,  the  procedure  preferred  here  is  to  recast 
the  problem  in  a form  such  that  the  two- indepen- 
dent-variable dependence  is  reduced  to  dependence 
on  a single  independent  variable,  time.  This 
may  be  achieved  by  use  of  the  method  of  weighted 
residuals,  or,  more  simply,  by  use  of  a straight- 
forward integral  method  [28];  the  result  is  an 
initial-value  problem,  consisting  of  simulta- 
neous quasilinear  ordinary  differential  equations 
in  time,  with  appropriate  starting  conditions. 

The  adequacy  of  this  approximate  procedure 
depends  upon  the  ability  to  anticipate  the 
spatial  variation  of  the  dependent  variables. 

An  alternative  method  to  reduction  to  a one- 
independent-variable  problem  is  to  introduce  the 
integral  transform  over  space,  with  the  inversion 
leading  to  two  simultaneous  nonlinear  Volterra 
integral  equations  [26]  for  H(t),  f(t);  while 
this  reduction  is  formally  exact  for  cases  of 
interest,  convergence  of  the  iterative  technique 
required  for  solution  of  the  integral  equations 
must  be  achieved. 

For  completeness,  it  may  be  remarked  that 
the  simplification  of  the  Stefan-problem  formu- 
lation is  based  on  the  smallness  of  the  reaction- 
zone  thickness  relative  to  the  preheat-zone 
thickness,  for  large  dimensionless  Arrhenius 
activation  temperature.  Phenomena  in  the  burned- 
gas  region  are  taken  to  occur  on  the  diffusive 
scale  of  the  preheat  zone.  Now,  in  fact,  a 
third  spatial  scale  implicit  in  the  adoption  of 
the  range  has  been  introduced;  in  the 

Cartesian  coordinate,  reference  may  be  made  to 
the  dimension  of  the  "container"  L*;  if  the 
Peclet  number  Pe  based  on  adiabatic  flame  speed 
is  large,  i.e.,  if  Pe  = L*/(D*/u*)  » 1,  then 

even  the  diffusive-convective  preheat  and  post- 
flame domains  become  narrow  relative  to  scale  L* 
In  this  circumstance,  the  diffusive  terms  may  be 
dropped  from  the  Stefan-problem  formulation, 
provided  the  effect  of  the  flame  is  appropriately 
accounted  for  in  the  resulting  nondiffusive 
problem.  In  such  a nondiffusive  formulation,  the 
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dependent  variables  themselves  (except  for  the 
pressure)  are  discontinuous  across  the  parametri- 
cally-inserted flame.  Further  details  on  this 
level  of  approximation  are  deferred  to  later 
publications;  here,  results  for  the  already- 
developed  Stefan  problem,  retaining  the  diffusive 
scale,  are  obtained. 

PROBLEMS  IN  FLAME  PROPAGATION 

TRANSIENTS  - Before  problems  including 
flame/wall  interaction  are  undertaken,  some 
unsteady  flame  propagation  problems  in  the 
infinite  domain  -°°  <_ip  <_  « are  examined  to  eluci- 
date the  formulation,  the  methods  of  solution, 
and  some  properties  of  the  phenomena  under  study. 
Here  the  integral  method  is  used  to  study  the 
effect  of  Lewis-Semenov  number  on  reattainment 
of  steady  flame  propagation  after  a small  heat 
extraction  from,  or  addition  to,  the  flame. 

The  isobaric  initial/boundary-value  problem 
is  given  by  Eqs.  34-44  and  Eq.  22,  with  p = 1, 
ik  ■+  -°°,  ^ = 0.  The  boundary  conditions 

(Eqs.  45-47)  become 


h - 1 »s  * - - (for  t i 0);  H - 0.  V - I as  « - — (for  t l 0). 


(50) 


In  anticipation  of  an  integral  method  of 
solution,  results  are  sought  in  the  forms,  com- 
patible with  Eq.  39, 

h • 1 - (1-fi)  c ) for  f ' * < ■,  t ! Oi  (51 ) 


H ■ R exp(-(y-*)/le  b). 


(52) 


Y * 1 - exp(-(*-*)/d) 


for  -<*<»,  tiO,  (53) 


with  c,  b,  d,  as  well  as  H and  'F,  = fncs(t) 
These  forms  are  compatible  with  the  initial 
conditions  (Eqs.  42-44)  if  c -*■  Cg,  b -*■  bQ, 

d -*•  dg,  and  H -*■  fig  as  t ■*  0.  Insertion  of 

Eqs.  51-53  in  Eqs.  40-41  yields 


i 

o 


. ♦ i 


MJ.l 

c b d 


(54) 


At  t = 0,  for  given  values  of  Hg  and  Cg,  these 
relations  give  bQ  and  dQ;  for  H(t)  = 1,  c(t)  is, 
in  general,  irrelevant.  If  Hg  is  close  to  unity, 
so  are  bg  and  dg. 

The  evolution  in  time  of  the  solutions 
(Eqs.  51-53)  requires  the  finding  of  b(t),  c(t), 
d(t),  H(t),  and  Y(t).  Two  algebraic  eqquations 
are  obtained  by  the  substitution  of  Eqs.  51-53 
in  integrals  of  Eqs.  35-37.  For  example,  if 
prime  denotes  ordinary  derivative  with  respect 
to  time,  from  Eq.  36 


/K*-/i?‘-(fi)-  (55a) 


/ hi-  “*“!]  [?  I1’  - (»-*)d  jj  d*  ■ - g — > (55b) 


d(t*  - d ) ■ -1.  (55c) 

Similarly,  Eq.  35  gives 

H'  - 0 - R)  [tfc  4 r] - • (56) 

and  Eq.  37  gives 

“>’0r-n T"r)"  (57) 

If  Eq.  54  is  enforced  in  algebraic  form  at  t = 0 
to  determine  bg,  dQ,  then  the  partly  differential 

set  is  conveniently  rendered  a set  of  coupled 
quasilinear  ordinary  differential  equations  by 
differentiation  of  Eq.  54: 


♦ D ♦ 


o. 


- 0. 


(58) 

(59) 


For  Hg  = H = 1,  Eq.  56  is  an  identity,  c plays 
no  role,  d = dQ  = 1 from  Eq.  54  and/or  Eq.  58, 

b s d = 1 from  Eq.  54  and/or  Eq.  58,  and 
'F1  = -1  from  Eq.  55c  and/or  Eq.  57,  so  that 
V = -t,  since  'Fg  = 0.  Thus,  as  a special  case, 

this  set  can  recover  steady  laminar  flame  propa- 
gation (here  expressed  in  a laboratory,  as 
opposed  to  a flame-fixed,  frame  of  reference). 
Since  H = 1 for  appreciable  reaction,  Eq.  58 
gives  interesting  insight  into  the  relative 
importance  of  the  reaction  order  and  the  non- 

dimensional  activation  temperature  B.  Incremen- 
ting B by  two  has  the  same  effect  on  the  solution 
as  incrementing  Vp  by  unity  for  H = 1;  since 

B = 0(10),  Vp  = 0(1),  typically,  an  appreciable 

change  in  Vp,  such  as  an  increase  from  one  to  two, 

has  modest  effect  on  the  results. 

The  following  results  concerning  the  effect 
of  Lewis-Semenov  number  Le  and  of  Hg  on  the  rate 

of  recovery  of  laminar  propagation  speed  are 
obtained  by  numerical  integration  of  Eqs.  55c, 
56-59,  for  Cg  = 1,  3 = 10,  Vp  = 1.  For  Le  = 1, 
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Hq  * 0.9  (such  that  Tq  * -0.25)  results  in 

H * 0.97  and  ¥'  = 0.9  for  t = 10,  but  Hq  » 0.5 

(such  that  Tq  = -0.005)  results  in  effectively 

no  recovery  at  all,  even  by  t = 100.  Hence,  for 
too  large  a perturbation,  steady  flame  propaga- 
tion effectively  is  not  recovered,  according  to 
the  model.  Also,  for  Le  * 1,  Hq  * 0.99  (such 

that  Yq  = -0.93)  results  in  ¥'  = -0.99  for 

t * 3,  but  Hq  = 1.01  (such  that  Yq  = -1.07) 

results  in  = -1.01  for  t = 1.9;  hence, 
recovery  is  faster  for  flame-temperature  augmen- 
tation than  for  flame-temperature  depression. 

For  Hp  = 0.9,  at  Le  = 0.7  (such  that  Yq = -0.26), 

H = 0.99  at  t = 17.5  and  f = -0.99  at 
t * 23.5;  for  Hq  = 0.9,  at  Le  = 1.3  (such  that 

^ = -0.25),  H = 0.99  at  t = 7.9  and  i,'  = -0.99 

at  t = 6.5  Thus,  for  larger  Lewis-Semenov 
number,  both  the  flame  temperature  and  the  flame 
speed  are  recovered  sooner.  It  should  be  noted 
that  the  influence  of  the  Lewis-Semenov  number 
on  the  adiabatic  flame  speed  itself  is  implicit 
in  the  nondimensional ization;  the  influence  of 
the  Lewis-Semenov  number  on  the  rate  of  recovery 
of  steady  flame  propagation  after  an  initial 
temperature  perturbation  is  the  matter  under 
study. 

STRATIFICATION  - He"e,  the  Fourier  integral 
transform  is  used  to  derive  Vol terra  integral 
equations  describing  the  effect  of  nonuniform 
stoichiometry  of  the  premixture  on  laminar  flame 
propagation;  in  particular,  the  error  incurred 
by  using  locally  an  appropriate  steady-laminar- 
flame-propagation  formula  to  obtain  the  accelera- 
ting or  decelerating  flame  speed  is  sought.  Of 
course,  nonuniform  stoichiometry  of  the  fuel/air 
premixture  is  of  current  automotive  interest 
because  of  stratified-charge  engines. 

The  initial/boundary-value  problem  posed  by 
Eqs.  34-47,  and  Eq.  22,  is  recast  to  allow  a 
more  general  initial  condition,  and  to  incorpo- 
rate the  conditions  (Eqs.  40-41)  in  the  partial 
differential  equations  (35-37).  Thus,  over  the 
domain  <_  ip  <^®,  t ^ 0, 


>Y  ,% 

tt  - — J • -*  MV  • Y / , 

H 

(60) 

It  * Ly  —7  • f My  * 0, 

(61) 

where  <$(z)  is  the  Dirac  delta 
again  denotes  the  position  of 
definition. 

function,  'V(t) 
the  flame,  and,  by 

(62) 

where,  again,  ff  = H(t)  denotes  the  flame 
temperature. 

As  t -►  0, 


R 


Y 


T 


1.  f - Yft 


0.  i‘ 


1 • e*P(w)  ♦ A S(-v-t) 
0 


exp(*/le) 

I 


, (63) 

for  , ; o ( 64a ) 

for  v ; o.  (64b) 

for  * ; o. 

(65) 

for  » ’ O, 


where  S(z)  is  the  Heaviside  unit  steo  function. 
The  step-type  increment  (or  decrement)  to  the 
upstream  unburned  premixture  is  characterized  by 
a magnitude  A at  distance  L from  the  initial 
flame  site  in  the  streamfunction  coordinate 
The  larger  L,  the  more  the  step  function  diffuses 
to  a smoothed-out  transition  from  Y = 1 to 
Y « (1  + A)  before  the  flame  arrives. 

The  boundary  conditions  are 

* * <1  ♦*),  T - 0 as  «.  - — (for  t > 0);  y - 0,  T I IS  « - * (for  t - 0) . 

(66) 

Clearly  the  flame  propagates  from  its  initial 
site  at  the  origin  in  the  direction  of  negative 

The  Fourier  transform  pair  over  the  spatial 
coordinate  is  introduced: 


7(s,t)  ■ f v(*.t)  *xp(-s<)d,  y(n.t)  ■ / Y(s.t)  exp(su )ds,  (67) 


with  the  transform  parameter  s = if..  For  brevity, 
only  the  mass-fraction  equation  is  developed;  the 
temperature  equation  is  readily  seen  to  be 
similarly  treated,  so  only  the  result  is  stated. 

If  tentatively  the  function  F is  formally 
treated  as  known,  and  if  contributions  from  the 
integrations  by  parts  are  ignored  (because  the 
contour  of  the  inversion  integral  with  respect 
to  singularities  on  the  real  axis  is  to  be  so 
chosen  that  the  boundary  conditions  are  incor- 
porated), then  application  of  the  transform  gives 

||  - V • -F  *«P(-ST);  (68) 

standard  introduction  of  an  integrating  factor 
for  this  first-order  differential  equation  gives 


7(s,t)  * V(s.O)  *xp(s?t)  - r f(t')  »«p  jt^(t-t’)  - 5T(t  ij  dt' . (69) 


Further,  the  taking  of  the  transform  of  the 
initial  condition  gives 


mmmm 
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?(«.0)  " • J * JTT  • J •«P(l»).  (70) 

Inversion  of  the  contour  integrals  [29]  gives 


The  physical  role  of  each  term  is  now  noted. 

Each  of  the  first  three  terms  on  the  right-hand 
side  of  Eq.  71  gives  the  time  evolution  of  one  of 
the  three  terms  on  the  right-hand  side  of  the 
initial  condition  (Eq.  64a);  the  first  term  of 
Eq.  71  is  derived  from  the  first  term  of  Eq.  64a, 
etc.  The  fourth  term  on  the  right-hand  side  of 
Eq.  71  is  consumption  at  the  thin  flame. 

The  result  of  Eq.  71  is  given  in  laboratory 
coordinates;  it  is  convenient  to  transform  to 
flame-fixed  coordinates: 


(*,t)  - (;.t).  with  c • * ♦ m(t). 


(72) 


Eq.  71  is  rewritten: 


The  flame  is  situated  at  c,  = 0 for  all  time  t. 
Enforcement  of  Eq.  39  gives 


0 • J *r,c(-  - 7 j'r,c(?7iT7^)i  (74) 


For  A = 0,  F = 1 , m = t is  a solution  of  this 
integral  equation;  for  A / 0,  H/l.soFjM, 
and  H and  m must  be  found.  Whereas  Eq.  74  gives 
one  equation  for  the  two  unknowns,  the  second  is 
obtained  from  consideration  of  Eqs.  61,  62,  65 
and  66.  It  is  found  that 


Enforcement  of_Eq.  39  gives,  if  it  is  recalled 
that  h(i|/,t)  -+  H(t)  as  tp  -*■  H'(t), 


For  H = 1 so  that  F=l,m=tisa  solution  of 
this  integral  equation  for  all  finite  values  of 
Le.  Eqs.  74  and  76  constitute  the  coupled 


relations  for  H and  m;  once  these  two  unknown 
functions  of  time  are  determined,  substitution  in 
Eqs.  73  and  75,  or  then-straightforward  numerical 
integration  of  Eqs.  60-66,  yields  the  mass- 
fraction  and  temperature  profiles  at  any  U,t). 

Because  a convergent  iterative  procedure 
for  the  treatment  of  Eqs.  74  and  76  is  not  an 
entirely  trivial  matter,  some  discussion  is  given 
here,  prior  to  presentation  of  numerically- 
obtained  results.  First,  the  integral  term  in 
Eq.  76,  denoted  I2[t,m(t);Le,vc,B] • and  the 

integral  term  in  Eq.  74,  denoted  Ij[t,m(t); 

Le,vF,B]  - I2[t,m(t);l,vp,8]  are  integrated  by 

use  of  linear  interpolation,  as  now  indicated. 
Discrete  time  intervals  of  interest  are  selected: 
ts>  s = 1,2,3,  ....  where  t$+j  > t$  and  tj  - 0. 

Then,  for  example. 


(77) 


V 

'j 

where  tl  is  a dummy  variable  of  integration.  Here 

J 


i [»(t  J -i»(t:)]  j 

Qjtxj  .H(  t j ) ."•{  * j ) ;ts.">(  ) ;Le,vp.e]  = F[tj,H(tj);\y,f]  exp  j - TTeTt~H’T( 

i V’  ( s(,'H(ti,)l  ( Wt  ) - »(t;n2 » 

' 3 I 2 Hltp  ( I 1 L*(,s  V ' 


; 9j[tj.ff(tj),«!(tj);ts.ni(ts);le,\>r,e]  (70) 

. [(t . , > ♦ (t ..  t;!i 

- ij  j s j 


over  the  domain  t.  < tl  < t.j 


Substitution 


M - Cj  - j+r 

of  Eq.  78  in  Eq.  77  and  integration  give 
I2[ts jLe.Vp.B]  as  an  algebraic,  as  opposed 


to  integral,  expression.  For  Le  = 1,  Eqs.  74 
and  76  sum  to 


H(ts)  ■ 1 ♦ 


(79) 


Use  of  Eq.  79  in  Eq.  76,  together  with  the  con- 
cepts of  Eqs.  77  and  78,  yields  one  transcen- 
dentally  nonlinear  algebraic  equation  for  one 
unknown,  m(ts).  A first  guess  may  be  obtained 

by  linear  extrapolation  from  m(t$  j)  and  m(tg  2); 

for  early  times,  m(t)  = t furnishes  guidance.  A 
second  guess  may  be  obtained  by  perturbing  the 
first  guess  by  a small  amount.  Newton's  rule, 
with  finite  difference  approximation  of  the  first 
derivative  based  on  the  last  two  iterates  from 
m(ts),  quickly  yields  effectively  invariant  suc- 

cesive  iterates  for  m(t$)  — and,  by  use  of  Eq.  79, 
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*or  H(t$)  --  and,  hence,  quickly  achieves  con- 

vergence to  a solution  for  t&.  Now,  let  s-*-(s+l), 

i.e.,  increment  time,  and  repeat  the  procedure. 
For  Le  t 1,  an  analogue  to  Eq.  79  is  obtained  as 

follows:  multiply  Eq.  74  by  (Le)'^2  and,  then, 
subtract  Eq.  76  from  Eq.  74  to  produce 


♦ {lz[ts.m(ts);L«,vr,6]  - "T77  I?[ts>n<ts,il*vF’6j} 


(80) 


The  (small)  difference  of  the  last  line  of  Eq.  80 
is  evaluated  as  a number,  by  use  of  the  best 
current  estimates  for  m(ts),  H ( ts ) ; the  first 

three  lines  of  Eq.  80  are  left  as  written.  Then 
Eq.  80,  for  Le  f 1,  serves  the  analogous  role 
that  Eq.  79  serves  for  Le  = 1,  i.e.,  substitution 
of  Eq.  80  in  Eq.  76  gives  one  nonlinear  algebraic 
equation  for  m(ts),  to  be  solved  by  iteration 

(using  Newton's  rule  to  obtain  the  root)  until 
invariance  of  successive  iterates  for  m(ts) 

indicates  convergence.  Use  of  this  result  for 
m(ts)  in  the  first  three  lines  of  Eq.  80  yields 

an  updated  H(ts).  These  values  for  m(ts),  H(ts) 

become  the  most  current  estimates,  and  the 
entire  procedure  can  be  repeated  until  conver- 
gence is  obtained. 

Results  are  presented  in  Figs. 2 through 
4;  spatial  profiles  for  Y and  T,  at  fixed  time 
t,  may  be  more  readily  obtained  from  numerical 

integration  of  Eqs.  60-66,  with  H(t)  and  y(t) 
known,  than  from  Eqs.  73  and  75.  A question 
arises  concerning  how  adequate  Eqs.  21  and  24 
are  as  a means  of  obtaining  the  instantaneous 
flame  speed  when  the  flame  is  propagating  through 
the  stratification  gradient,  such  that  the  effec- 
tive upwind  mass  fraction  varies  from  unity  at 
t = 0 to  (1  + A)  at  long  time.  The  following 
calculations  indicate  that  the  "quasisteady" 
calculation  Is  but  about  5%  or  less  in  error, 
and,  thus,  fully  adequate  for  most  engineering 
applications.  For  isobaric  flow  with  $ « 1, 

since  0 ^ ^Fu^"1,  ^rorn  2* 


0 5.0  10.0  15.0  20.0 


t 


Fig.  2a  - From  Volterra-intearal-equation  for- 
mulation of  thin  flame  propagation  through  a 
stratification  in  stoichiometric  ratio  of  the 
unburned  premixture,  the  dimensionless  flame 
temperature  H as  a function  of  time  t,  for  values 
of  the  dimensionless  activation  temperature  6 
and  stoichiometric  coefficient  for  fuel  Vp. 

Here,  Le  = 1,  A = 0.2,  L = 10.  It  may  be  noted 
that  R(0)  = 1 


(vf-1)/2  (vF+l )/? 

m-m  ft)  -ft (•?)! 


since  (e^Bj)  = (Ypui/Ypu2)-  Subscript  1 denotes 

the  initial  condition  (for  which  the  upwind 
asymptotic  value  of  the  mass  fraction  Ypuj  = 1); 

subscript  2 denotes  the  effective  upwind  asymp- 
totic value  of  Ypu  at  some  later  time.  Since 

m'(t)  is  a dimensionless  flame  speed  relative  to 
the  gas,  for  fixed  cold-premixture  density. 
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' 


where  mj  « 1.  The  largest  (or  smallest)  value 

that  the  ratio  (Ypu2/YFul)  achieves  Is  (1  + A), 

a value  Independent  of  both  L and  le;  this 
limiting  value  occurs  for  t + «>: 

■j(t  • •)  » (i  « a)  f •«[> • (83) 

Some  appreciation  for  the  relative  roles  of  para- 
meters Is  obtained  from  the  following  calcula- 
tions (unless  otherwise  specified,  p ■ 10, 
vF  ■ 1,  A • 0.2  so  m2(t  -►  «>)  ♦ 2.76):  for  A ■ -0.2, 

m^( t ► *•)  ♦ 0.23;  for  p * 15,  m^t  •*  »)  ► 4.19; 

for  vF  * 2,  m2(t  »>)  •*  3.31. 

For  Le  ■ 1.  an  Integral  exists  for  steady 
I laminar  Isobarlc  adiabatic  flame  propagation  In 
the  form  (Y  ♦ T)  » const.,  where  T ■ 0 at  the 
cold  boundary  and  Y ■ 0 at  the  hot  boundary. 
Heurlstlcal ly,  this  relation  can  be  employed  to 
define  an  Instantaneously  equivalent  steady 
problem  for  a flame  propagating  through  a fuel- 
mass-fractlon  gradient.  Specifically,  fuel  mass 
fraction  upwind,  where  the  normalized  temperature 
Is  zero.  Is  associated  with  the  Instantaneous 


temperature  at  the  flame,  where  the  fuel  mass 
fraction  vanishes:  YFu2(t)  i H(t).  For  t -*  0, 

YFu2  + ‘ for  t “•  yfu2  *■  ^ £ + A)- 

With  FFu2  = Ff(t),  a "quasisteady"  m^  may  be 

computed  from  Eq.  82c  to  compare  with  the  Instan- 
taneous m2  calculated  from  the  two  simultaneous 

Integral  equations.  For  the  parameter  ranges 

P ■ 0(10),  vF  ■ 0(1),  A ■ OdO'1),  the  flame 

speed  m'  Is  underestimated  by  no  more  than  61 
by  the  "quasi steady"  approach. 

For  Le  ■ 0(1),  the  equivalency  Ffu2(t)  H(t) 

may  be  defined  again  for  purposes  of  comparing 
a "quasisteady"  flame  speed  with  the  Instanta- 
neous mi  calculated  from  solution  of  a flame  . 

propagating  through  a premixture  of  nonuniform 
stoichiometry.  The  cited  Integral,  (V  ♦ T)  ■ 
const.,  no  longer  holds,  so  that  the  concept 
that.  In  steady  laminar  adiabatic  Isobarlc  flame 
propagation,  the  hot-boundary  temperature 
achieved  from  a given  oold-boundary  premixture 
Is  Independent  of  Lewis -Semenov  number.  Is 
employed.  However,  as  Le  departs  from  unity, 
the  accuracy  degenerates. 


Fig.  3a  - From  Volterra-lntegral-equatlon  for- 
mulation of  thln-flame  propagation  through  a 
stratification  In  stoichiometric  ratio  of  the 
unburned  premixture,  the  dimensionless  flame 
temperature  H as  a function  of  time  t,  for 
several  values  of  the  Lewls-Semenov  number  Le. 
Here,  A ■ 0.2,  L ■ 10,  vF  ■ 1,  p ■ 10.  It  may 
be  noted  that  H(0)  ■ 1 
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Fig.  3b  - From  Vol terra-integral -equation  formu- 
lation of  thin-flame  propagation  through  a strat- 
ification in  stoichiometric  ratio  of  the  unburned 
premixture,  the  dimensionless  time  rate  of  change 
of  the  flame  position  in  von  Flises  coordinates 
dm/dt  as  a function  of  time  t,  for  several  values 
of  the  Lewis-Semenov  number  Le.  Here,  A = 0.2, 

L = 10,  vp  = 1,  6 = 10.  It  may  be  noted  that 
m(t)  = -'t',  m(0)  = 0,  and  dm(0)/dt  = 1 


PROBLEMS  IN  FLAME/WALL  INTERACTION 

ISOBARIC  COLD-WALL  QUENCHING  - The  isobaric 
treatment  of  que  ching  of  combustion  near  a cold 
wall  is  more  pertinent  to  automotive  applications 
than  first  appears,  in  that  there  is  no  signi- 
ficant pressure  rise  owing  to  compression  over 
the  time  interval  for  flame/wall  interaction, 
for  not-too-lean  premixtures.  There  may  be 
significant  pressure  transients  owing  to  im- 
pressed piston  motion.  However,  for  lower  rates 
of  engine  revolution,  the  present  case  should  be 
adequate. 

For  the  domain  of  interest  0 <_  <1/  < °°,  t > 0, 
the  initial /boundary  value  problem  becomes 
Eq.  22,  with  p = 1,  together  with 

Y = 0.  f£  - U - 0 for  r < * < - ; (®4) 

1 Hi 

(85) 


At  the  cold  noncatalytlc  wall,  the  boundary  con- 
ditions are 


Fig.  4a  - From  Vol terra- integral -equation  formula 
tion  of  thin-flame  propagation  through  a stratifi 
cation  in  stoichiometric  ratio  of  the  unburned 
£remixture,  the  dimensionless  flame  temperature 
H as  a function  of  time  t,  for  values  of  the  mag- 
nitude of  change  in  unburned  fuel  mass  fraction  A 
and  the  dimensionless  initial  distance  of  this 
change  from  the  flame  L.  Here,  Le  = 1,  vp  = 1, 

B = 10,  11(0)  = 1.  It  may  be  noted  that  the  case 
A = -0.2,  L = 10  has  been  terminated  at  a time 
beyond  which  there  is  little  of  interest 


H - o.  !j  - 0 <s  « - 0 (for  t : 0).  (86) 

while  at  the  hot  reservoir. 


► * 1 «S  * * ■ Ifor  | ; 0)  187) 


At  the  flame,  = f(t),  given  initially  by 
V(0)  = Vq  » I (given  const.),  conditions  of  Us 

40  and  41,  with  p * 1,  hold;  the  quantity  F ( t , 
introduced  in  Eq,  62,  is  to  be  used  again  below 
The  initial  conditions  are  given  by  Eqs.  43  and 
44,  with  Hq  * 1,  dp  * 1,  bg  * 1.  A fuel -lean 

thin  vigorous  flame,  initially  far  from  the  cold 
impervious  noncatalytlc  wall,  propagates  towards 
the  wall  and  interacts  with  It. 

Derivation  of  Vol terra  ntegral  Equations 
by  Integral  Transform  - FTr  ‘ ~ the  integral - 
equation  reduction  to  a problem  in  time  only  is 
executed.  The  semi-infinite  domain  in  is 
extended  to  the  full-infinite  domain  via  odd 
symmetry  for  the  temperature,  and  even  symnetry 
for  the  fuel  mass  fraction,  about  ip  * 0;  the 
boundary  conditions  at  the  thin  flane  4/  * f(t) 


- ..  - ^ 
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are  incorporated  In  the  partial  differential 
equations  by  use  of  the  Dirac-delta  function,  as 
was  done  earlier.  Hence, 


- F4[*  - 1)  - f«c*  ♦ »].  (88) 

It  - Ip  • Fit*  - »]  - '«[*  4 »].  (89) 

31  3* 

Application  of  the  Fourier  integral  transform 
over  i *>  is  now  straightforward.  However,  the 
initial  conditions  must  be  modified  slightly  to 
preserve  the  even  symmetry  of  Y(vp,t) , and  the  odd 
synmetry  of  T(ip,t).  Here,  it  is  taken  that,  as 
t 4 0, 


Y 

0 


0 


cosh(a  Yq)  - coshfa  *) 
cosMj  »0)  - 1 


For  |»|  > ♦„ 


for  |*|  4 Y0. 


(90) 


for  |»|  > »0 


for  |»|  . »0> 


(91) 


where  a, 6 are  constants  to  be  assigned,  such 
that  these  forms  simulate  initial  conditions  at 
the  flame  front  (i|<  ♦ Yq  ) as  closely  as  possible. 


Satisfaction  of  the  boundary  conditions  for 
i*>  * vQ  for  an  adiabatic  flame  gives 


dY„  a Sl"M»  »„) 

ar  ’ • ’ 1 

(92) 

» s2  - ' — » «•»"<»  V ■ a ■ 

(93) 

(94) 

Also,  TqU)  and  YQU)  are  virtually  identical  to 

the  initial  profiles  first  discussed. 

As  before,  the  Fourier  transform  of  Eq.  88 
over  + yields  formally 


3 r (95) 

*(*.«)  • *(».01  np(,;i) 

I 

/ f(l  > j«»p[-iY(t')]  4 P»pfi*(«')]}  dt'. 

(96) 

But,  from  Eq.  90,  with  A ■ cosh(nT0)/[cosh(aS'0)-l] , 


YU.O)  ■ « 


| ""Mi  *„)  *lnnf (»-a)  »0) 
j l ' ffs-a)  coih(n  ♦„) 


j(nh[(s*a)  »0)  j 
fT*4")  COlhfa  »0)  | 


• (97) 


Insertion  of  Eq.  97  in  Eq.  96  and  inverting  gives 


Fig.  4b  - From  Vol terra-integral  equation  formula- 
tion of  thin-flame  propagation  through  a stratifi- 
cation in  stoichiometric  ratio  of  the  unburned 
premixture,  the  dimensionless  time  rate  of  change 
of  the  flame  position  in  von  Mises  coordinates 
•dm/dt  as  a function  of  time  t,  for  values  A and 
L.  Here,  Le  = 1 , vp  = 1 , 8 = 10,  m(0)  = 0, 
dm(0)/dt  =1.  It  may  be  noted  that  the  case 
A = -0.2,  L = 10  has  been  terminated  at  a time 
beyond  which  there  is  little  of  interest 


A(K'r,&)',rf(^)] 


i r 

- 1 CoYhTal0l 

■'O'  • "'") 

i r 

• 1 eo*KM07  [»' 

’ """""o'  i v?f'‘  / \?t”‘  f\l 


(98) 


The  terms  which  multiply  A describe  the  time 
evolution  of  the  initial  conditions,  while  the 
integral  term  describes  chemical  consumption. 
Similarly, 


_ 


The  conditions  at  the  flame,  o = 1,  are: 





1 #Lt  6 t ♦ 6* 
’ I »<nh(B  fQ) 


W— 

L \2(L« 


TT77  * 8(L*  t)V? 


?)'*r,(S^tB(utl'/2)] 


h - R.  H - R.  Y - 0. 


(104) 


• i i'w'v1  [•"  (jtifa  ■ ■ •r,(sri?7! ' ““  [(g . L,  a)  - u g]  g . ' 

tgSflL..  m 

J("l«>12-{  [ (t  - f)1/2  V ' 


6 1-R  ( 

■?Tr 


• 10-1 


(for  t > 0). 


(105) 


Enforcement  of  Eq.  39  gives  two  simultaneous 
nonlinear  Vol terra  integral  equations  for  the 
two  unknowns  ^(t),  H(t)  from  Eqs.  98,  99;  once 
solution  is  obtained,  spatial  profiles  may  be 
found  from  the  same  two  equations.  Of  course, 
by  calculation  of  normalized  residuals  for 

T[¥(t),t][«H(t)]  and  Y[¥(t),t],  Eqs.  98  and  99 
also  may  be  used  to  characterize  the  accuracy  of 
solution  for  ff(t)  and  yft)  obtained  by  an  alter- 
native, approximate  method,  such  as  the  integral 
method  now  discussed. 

Derivation  of  Initial -Value  Problem  by 
Integral  technique  - the  problem  is  given  by 
partial  differential  equations  (Eqs.  34-85); 
boundary  conditions  (eqs.  86-87);  supplementary 
relations  (Eq.  22)  and  flame  conditions  (Eqs. 
40-41),  all  with  pjf  1;  and  initial  conditions 
(Eqs.  43-44,  with  Hq  = 1,  dg  = 1,  bg  = 1). 

First,  the  (so-called)  Landau  coordinate 
transformation  [30]  is  introduced  to  fix  the 
flame  position: 


(*.t)  - (o.t),  with  o ■ */Y(t).  (100) 


The  initial  conditions  are,  for  H(0);Hq=1, 

f(0)  =Yg  » 1: 

h - hQ  • 1 as  t ->  0 (for  1 < o < -);  (106) 


H - Hq  ■ exp(-Y0(l-o)/Le),  Y - Yq  • [1  - exp {-fQ( 1 -o) >]  »s  t - 0 

(for  0 i o i 1). 

(107) 

For  4*q  » 1,  as  noted  previously,  these  condi - 
tionsuare  compatible  with  the  boundary  conditions 
to  within  (exponentially)  small  error. 

The  following  representations  of  the  depen- 
dent variables  are  taken  to  be  good  approximations 
for  all  time: 


h(o.t)  ■ 1 - (i  - R(t))  »p  j-  3^11$! | for  i : o : ■ ; (108) 


The  initial/boundary-value  problem  for  Y(o,t), 

h(o,t),  H(o,t)  is  restated  (here  primes  denote  # t n , s1nh { ^ 

ordinary  derivatives  with  respect  to  time,  H(g,t)  * ff(t)  ~<  Y(tr~r  (109) 

denoted  by  t):  UTBmr 


Y = 0, 


ah  r 

3t  ‘ Y 


ah  Le  a2h 

* ‘ 7 1? 


0 


for  !<<?<•; 


(101) 


Y(o,t)  - 1 


cosh  { af 

»} 

f or  0 : o ; \ . 


(110) 


a y r „ aY  in, 
at " t 0 ao  ' It  7T 


n ah  y*  an 

°-  at  ' r°  a5  ' 


’ ° for  0 * 0 « 1 (102) 


The  boundary  conditions,  for  a cold  Isothermal 
noncatalytic  wall,  and  for  the  far  field,  are 


For  b(0)  = bg  = 1,  d ( 0)  = dg  = 1,  these  adopted 

forms  are  closely  compatible  with  the  initial 
conditions.  The  two  flame  conditions  at  o = 1 
involving  gradients,  and  the  three  partial 
differential  equations,  yield  five  relations 
for  the  five_unknown  functions  of  time  b(t), 
c(t),  d(t),  H(t),  v(t).  Initial  values  are: 


0 11  » • 0 (for  t : 0);  h < I »!  O » ■ (for  t i 0). 


(103) 


R(0)  = Hq  ■ 1,  Y(0)  5 >>  1,  given  const.;  c (0)  5 cQ  • 1 (formally).  (ill) 


The  values  of  b(0)  bg,  d(0)  dg  are  to  be 

In  the  manner  now  described. 

The  flame  conditions  (Eq.  105)  give 


• • (WWn  •¥>- 


0 * ' ■ • (is)f  ‘ • (a)  u"*,T/d)>  (119) 


* • • (S)f  ■ » «J 'mi  n • (£).  • h M&i 


t««n(jr/a)  . Rvr*’  t«p  j-  | IHlJ 


At  t * 0,  for  fg  » 1,  these  give  bg  dg  » 1 

by  inspection.  Differentiation  with  respect  to 
time  gives  for  Eqs.  112  and  113,  respectively. 


[H  I ul^/U  b)  ' ewl/l,  b)|]  7 * °‘R>  tS 

♦ t*nh(i/d)  . _ J&41-  la;  - f iial  . M 4t Pi 1 XI 

co»IT(»/d)  J d‘  [d  cosMY/d)  b »inlr(Y/le  b)J  ' 

4 H {bisiKimrsT 4 c {]  * 01  (ll 


r-  • • (£)f  ■ c) 


Hence,  from  Eqs.  116  and  1 1 7 


[(»/ d)  »«chJ(*/d)- tinb(i/d)]d‘  * (tti**(t/d)]v  • - tfB^t/dJ  , (122) 


1 ’ s'inhlwii  VT)]  **  ’ Le  tinh(«/le  k)]’  ‘4BR,*rft‘  (123) 


and  from  Eq.  118 


(l.R)c'  ♦ i»Jl  »■  . c r.  . lA 


' * rfnKft'hi  cotS(tM)  j d * | l(nbf!/jy‘coib7*7(f)  j r 

*>4,,4i]f ,#-  <115) 

In  general,  the  partial  differential  equa- 
tions may  be  integrated  with  respect  to  a,  upon 
imposition  of  the  boundary  conditions  and  of  the 
flame  conditions,  to  yield 


\ 

wh*r*  0 • f Y d.'.  fcn(t).  (\.t).  fcntth 


i 

Wh#r#  * " f M fnc^*  (|%)  * S ( 1 • t ) • fnc(t).  |:JH)  • (O.t).  fcn(t); 

■»  ' * (117) 

■v  (&),•«»-»  it. 

wbtr*  n • f (1-Mdo,  (£)  • i®  (l,»|,  (118) 

"1  * 

For  the  specific  forms  adopted. 


This  last  equation  states  that,  if  Hg  * 1,  H * 1 

at  all  subsequent  time;  this  artifact  of  the 
approximate  procedure  is  circumvented  by  starting 
with  a value  of  R slightly  displaced  from  unity. 
The  fact  that  the  calculational  procedure  tends 
to  establish  a steady  propagating  flame  while 
still  away  from  the  wall  (4’  » 1),  provided  Hg 

is  not  too  far  from  unity,  is  established  in  the 
section  on  the  effects  of  transients  on  flame 
propagation. 


Equations  114-115,  122-124,  subject  to  the 
Initial  conditions  implicit  in  Eqs.  111-113, 
constitute  the  sought-after  initial-value  problem 
for  b(t),  c(t),  d(t),  R(5),  H'(t).  Substitution 
of  the  results  in  Eqs.  108-110,  22  gives  the 
solution  for  Y(<ji,t),  T ( , t ) , H'(t);  results  depend 
on  the  parameters  K,  8,  Vp,  and  Le.  (tg,  cQ,  Hg 

must  also  be  selected,  but  there  is  minor  sen- 
sitivity to  these  quantities  introduced  by  the 
method  of  solution.)  (See  Figs.  5-8.) 

Other  output  quantities  of  interest  are  the 
rate  of  fuel  consumption  at  the  flame  (i.e.,  the 
derivative  of  the  fuel  mass  fraction  at  the 
flame),  the  fuel  mass  fraction  at  the  wall,  and 
the  temperature  gradient  at  the  wall;  these 
quantities  are,  respectively  (if  Xf(t)  is  the 

site  in  the  Cartesian  coordinate  of  the  flame 
t(t)): 
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(!o)f  ’ '(t/d>  UWiCr/O  ' '(,/d)  tanh(y/d)  — 


Fig.  5a  - From  integral -method  results  for  a cold 
wall  (at  temperature  T = 0),  profiles  for  the 
mass  fraction  Y as  a function  of  the  coordinate 
o[=  liV'f'(t)]  at  time  t = 8.04,  for  several  values 
of  the  dimensionless  Arrhenius  activation  tem- 
perature 8.  Here,  K = 6,  Le  = 1,  vF  = l,  and  the 
initial  time  is  defined  by  ^(t  = 0)  = fg  = io. 
This  is  the  first  of  five  figures  which  describe 
the  augmented  quenching  for  flame  interaction 
with  a cold  wall  with  increase  of  the  activation 
temperature  3 


Fig.  5b  - From  integral -method  results  for  a cold 
wall  (at  temperature  T = 0),  profiles  for  the 
temperature  T as  a function  of  the  coordinate 
o[*  i|/f(t)]  at  time  t = 8.04,  for  several  values 
of  the  dimensionless  Arrhenius  activation  tem- 
perature 8.  Here,  K = 6,  Le  = 1,  vF  = 1,  and  the 
initial  time  is  defined  by  '♦'(t  * 0)  ■ Hq  = 10 

I 


■ t - d Unh(t/d)  , 
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where  the  subtracted  term  In  Eq.  132  Is  Identi- 
fied to  be  the  depletion  of  fuel  ahead  of  the 
flame  owing  to  the  presence  of  the  flame. 

ISOBAR  I C ACCELERATION  TOWARD  AN  ADIABATIC 
WALL  - Derivation  of  Volterra  Integral  Equations 
by  Integral  Transform  - The  entire  derTvatYon 
leading  to  Fq.  WTiolds  without  modification. 
However,  because  Eq.  86  is  altered  to 


0 44  » * ® (*«r  1 


(133) 


Eq.  81  is  revised  to 

U ■ i«  LJ  • r »(*■»)  < r t ; o,  (134) 


to  enforce  the  evenness  of  T(d’ . t)  about  Y ■ 0 
for  the  adiabatic  case.  Also,  Eq.  92  is  revised, 
so  that,  as  t -»  0 


T(*.t)  - •„(*)  • 


\ 


•or  » 10 
•or  |»|  » •„ 


(135) 


Fig.  5c  - From  Integral -method  results  for  a cold 
wall  (at  temperature  T - 0),  profiles  for  the 
position  of  the  flame  In  von  Mlses  (mass-weighted) 
coordinates  Y as  a function  of  time  t,  for 
several  values  of  the  dimensionless  Arrhenius 
activation  temperature  8.  Here,  K ■ 6,  Le  - 1, 

F ■ 1.  and  the  Initial  time  is  defined  by 
Y(t  ■ 0)  • Y„  - 10 


Fig.  5d  - From  integral -method  results  for  a 
cold  wall  (at  temperature  T * 0),  profiles  for 
the  net  remaining  fuel  N (defined  by  Eq.  132) 
as  a function  of  time  t,  for  several  values  of 
the  dimensionless  activation  temperature  8. 
Here,  K = 6,  Le  = 1,  vp  * 1,  and  the  Initial 
time  is  defined  by  Y(t  * 0)  * Yq  * 10 


Fig.  5e  - From  Integral-method  results  for  a cold 
wall  (at  temperature  T « 0),  profiles  for  the 
flame  temperature  R as  a function  of  time  t,  for 
several  values  of  the  dimensionless  Arrhenius 
activation  temperature  8.  Here,  K - 6,  Le  - 1, 

VF  * 1.  and  the  Initial  time  is  defined  by 
♦(t  - 0)  - Y0  > 10 


Fig.  6a  - From  integral -method  results  for  a 
cold  wall  (at  temperature  T = 0),  profiles  for 
the  position  of  the  flame  in  von  Mises  (mass- 
weighted)  coordinates  t as  a function  of  time  t, 
for  several  values  of  the  Lewis-Semenov  number 
Le.  Here,  & = 10,  K = 6,  vp  = 1,  and  the  initial 
time  is  defined  by  t(t  = 0)  = Vp  = 10.  This  is 
the  first  of  three  figures  which  describe  the 
augmented  quenching  for  flame  interaction  with  a 
cold  wall  with  increase  of  the  Lewis-Semenov 
number  Le 


Fig.  6b  - From  integral-method  results  for  a 
cold  wall  (at  temperature  T = 0),  profiles  of 
the  net  remaining  fuel  N (defined  by  Eq.  132) 
as  a function  of  time  t,  for  several  values  of 
the  Lewis-Semenov  number  Le.  Here,  6 = 10, 

K = 6,  Vp  = 1,  and  the  initial  time  is  defined 
by  y(t  = 0)  = Tq  = 10 


Eq.  93  is  revised,  so  that 


s1nh(BYn) 


For  fg  » 1,  6 = (1/Le),  as  before.  Fourier 

transformation  over  ip , integration  over  time 
and  inversion  of  the  transform  gives 


coihfB't.) 


Fig.  6c  - From  integral -method  results  for  a 
cold  wall  (at  temperature  T = 0),  profiles  of  the 
thin-flame  temperature  R as  a function  of  time  t, 
for  several  values  of  the  Lewis-Semenov  number  Le 
Here,  6 = 10,  K = 6,  vp  = i,  and  the  initial  time 
is  defined  by  >f(t  * 0)  = * 10 
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Enforcing  T[Y(t),t]  = ff(t)  yields  an  integral 
equation  to  complement  the  one  found  by  enfor- 
cing Y[Y(t),t]  3 0 for  Eq.  98. 

Derivation  of  Initial-Va lue  Probl em_by 
Integral  JecTTnlque  - "The  entire  discussion  from 
Eq's.  10ff-Tl0¥oTds  unaltered,  except  for  the  y _ 

temperature  boundary  condition  at  the  wall;  i.e.,  * 

the  first  of  Eq.  103  becomes  133,  so  Eq.  109  is 
revised  to 

H(o.t)  • H(t)  c°‘h  ( 1381 


cosh  | 

COib  { iSt 


for  o : C i 1 . 


Thus,  Eq.  112  becomes 


"h^f/d)  t B unh(Vte  b)  _ (l-B)  , „ 


so  Eq.  114  becomes 


(T/l»  b)  s»ch?(»/Lf  b)j 


|h  | tanh( t/Le  b)  ♦ (r/l»  b)  s*ch?(*/l»  b)|j  b,  ♦ (l-B)  ^ 

|t»n*i(*/0)  * (t/d)  s»ch?(»/d)J  5,  . . t,  teJSA.HL  t*ch.' 

[B|j*aipilo!j]El.o.  (1 


♦ tinh(f/t»  b) 


Also,  Eq.  123  must  be  replaced  by 


r. (’t/it  i 

y * »1nh(t/Le  t)  < 


I t*nh(v/L»  b) 
If 


After  Eqs.  139  and  113  are  used  to  assign  bQ,  dQ 

for  specified  values  of  Le,  0,  Vp,  Hq,  and  Yq, 

then  straightforward  integration  in  time  follows 
by  use  of  Eqs.  115,  140,  122,  141,  and  124. 

Output  quantities  of  interest  for  the  adia- 
batic wall  are  the  rate  of  fuel  consumption  at 
the  flame  and  the  fuel  mass  fraction  at  the  wall; 
these  are  given  by  Eqs.  125  and  126,  respectively. 
The  temperature  at  the  wall  is  also  of  interest: 


• R secM*/L«  b).  (142/ 

Since  Y goes  from  yq  ( »1)  to  0,  Hw( t ) goes  from 

effectively  zero  to  the  temperature  at  the  flame 
as  the  flame  front  accelerates  towards  the  wall. 

The  association  of  Cartesian  and  von  Mises 
coordinates  is 


'Vbvcr  ,4 

/VV 

v ' \ /♦ 

V \ 
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Fig.  7 - From  integral -method  results  for  a cold 
wall  (at  temperature  T = 0),  profiles  for  the 
(normalized,  stoichiometrical ly  adjusted)  mass 
fraction  Y and  for  the  dimensionless  temperature 
T as  a function  of  the  dimensionless  Cartesian 
coordinate  x,  for  several  values  of  dimensionless 
time  t.  Here,  R = 10,  K = 6,  Le  = 1,  vp  = l. 

The  initial  time  corresponds  to  von  Mises  (mass- 
weighted)  spatial  coordinate  Y(t  = 0)  = Yg  = 10* 
Whereas  Y‘(0)  = -1,  the  initial  flame  temperature 
H(0)  = 1,  and  the  initial  net  remaining  fuel 
N ( 0 ) = 9,  the  calculation  is  terminated  at 
t = 8.99  with  Y‘  (0)  = 0.123,  H = 0.896,  and 
N = 0.434.  The  flame  speed  Uf  increases  from  an 
initial  value  of  -1.00  at  t = 0 to  a maximum  value 
of  -1.37  at  t = 7.42,  then  decreases  rapidly 
[with  uf(t  « 8.99)  = - 0.372 ] 


_ I ■(••i.xiraa;1! 
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where  Xf(t)  is  given  by 


Fig.  8 - From  integral -method  results  for  a cold 
wall  (at  temperature  T = 0),  the  speed  of  the 
flame  in  a laboratory  frame  of  reference  Up  as  a 
function  of  time,  for  several  values  of  the  exo- 
thermic! ty  parameter  K.  Here,  p = 10,  Le  = 1, 
vp  = 1,  and  the  initial  time  is  defined  by 
Y (t  = 0)  = Ya  = 10 


(w) 

uf  . ^ ■ (i  . « H)f  ♦ » H.  <B»mn  . (145) 

. li  WW'  »)  ; (146) 
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Results  are  given  in  Figs.  9-11. 

NONISOBARIC  ASPECTS  OF  FLAME/WALL  INTERACTION 

For  a contained  premixture,  the  isobaric 
approximation  requires  reconsideration,  and  the 
initial/boundary-value  problem,  as  presented  in 
Eqs.  22,  34-47,  is  to  be  examined.  The  integral  - 
equation  approach,  developed  above  for  isobaric 
interaction  of  a flame  with  both  a cold  iso- 
thermal wall  and  also  an  adiabatic  wall,  has 
been  generalized  to  the  nonisobaric  case  for 
both  types  of  wall  boundary  conditions.  The 
generalization  does  require  simultaneous  solu- 
tion of  six  quasil inear  coupled  ordinary  dif- 
ferential equations,  since  the  spatially 
invariant  pressure,  p(t),  joins  the  five  unknown 
functions  treated  for  the  isobaric  case.  Inte- 
gration of  the  middle  equation  of  Eq.  22  over 
the  flow  domain,  enforcement  of  the  impervious 
condition  on  the  velocity  at  the  containing 
walls,  and  substitution  of  the  adopted  spatial 
variation  of  the  temperature  yield  the  sixth 


Fig.  9 - From  integral -method  results  for  an 
adiabatic  wall,  profiles  of  the  (normalized, 
stoichiometrically  adjusted)  fuel  mass  fraction 
Y as  a function  of  the  dimensionless  Cartesian 
coordinate  x,  at  several  times  t,  where  the 
initial  position  of  the  flame  in  von  Mises 
coordinates  H'qMO.  For  this  nominal  case, 

Le  = 1,  K = 5.  Since  the  nondimensional  tem- 
perature T = 1-Y,  the  wall  temperature  T(0,t) 
rises  rapidly  to  the  adiabatic  flame  temperature 


equation  to  complement  the  five  obtained  from 
Eqs.  35-37,  40-41.  Details  and  results  are 
reserved  for  another  publication. 

The  conjecture  here  is  that  significant 
flame/wall  interaction  occurs  over  a few  diffu- 
sion scales  only;  thus,  as  long  as  the  bulk-gas 
flame  speed  is  appreciable  (i.e.,  the  premixture 
is  not  so  fuel -lean  that  the  laminar  propagation 
speed  is  highly  retarded),  the  quenching  of  the 
flame  near  a cold  wall,  or  the  acceleration  of 
the  flame  near  an  adiabatic  wall,  is  isobaric. 
The  compression-enhanced  pressure  level  attained 
by  the  burning  of  the  entire  cohfined  mass  of 
gas  is  the  pressure  level  to  use  in  the  near- 
wall analyses.  The  final  pressure  attained  in 
a homogeneous  fuel -lean  premixture,  uniformly 
at  pressure  p*  at  the  initiation  of  burning 

in  a container  of  length  L*,  may  be  shown  bv 
thermodynamic  arguments  to  be  p^  = P*  + (y-1)* 

YFu  P*  Q*  L*.  in  the  absence  of  losses  (which 

may  be  substantial).  The  anticipation  is 
that  results  obtained  from  an  inteqral- 
equation  method  for  a confined  premixture  will 
confirm  the  conjecture  that  flame/wall  inter- 
action is  virtually  isobaric,  as  far  as  compres- 
sion from  confinement  is  concerned. 

However,  there  may  be  significant  pressure 
change  over  the  time  interval  of  flame/wall 


Fig.  10a -From  integral -method  results  for  an 
adiabatic  wall,  profiles  of  the  dimensionless 
temperature  T as  a function  of  o(=  i^/f,  where 
is  the  von  Mises  spatial  coordinate  and  t is  the 
flame  position),  at  time  t = 7.52,  where  = 

YU  = 0)  = 10,  for  several  values  of  the  Lewis- 
Semenov  number  Le.  Here  the  dimensionless 
Arrhenius  activation  temperature  8 = 10,  the 
stoichiometric  coefficient  for  fuel  vp  = i,  the 
exothermicity  factor  K = 6.  Augmenting  the 
thermal  diffusivity  preheats  the  unburned  fuel 
more  effectively 
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Fig.  10b  - From  integral -method  results  for  an 
adiabatic  wall,  the  flame  position  in  von  Ilises 
coordinates  y as  a function  of  time  t,  for 
several  values  of  the  Lewis-Senenov  number  Le. 
Here,  3=  10,  Vp  = 1,  K = 6.  The  larger  the 
Lewis  number,  tne  more  rapid  the  final  accelera- 
tion (from  an  initially  constant-speed  propaga- 
tion) of  the  thin  flame  into  the  wall,  owing  to 
interaction  with  the  wall;  because  of  large 
terminal  gradients,  the  calculation  is  not 
completed  except  for  Le  = 1,  though  y -*•  0 at 
finite  time  in  the  model 


interaction  from  piston  motion  during  operation 
at  higher  rates  of  revolution  [14].  Prescribed 
movement  of  a wall  in  time  presumably  suffices 
to  simulate  piston  motion  in  the  one-dimensional 
unsteady  reacting  flows  under  study.  It  is  worth 
emphasizing  that  It  is  piston  motion  (and,  hence, 
wall  motion)  that  is  prescribed  as  a function  of 
time;  the  wall  motion  must  be  translated  into  a 
variation  of  pressure  with  time. 

DISCUSSION  AND  SUMMARY 

The  following  properties  of  results,  ob- 
tained by  use  of  a thin-flame  model  to  describe 
isobaric  unsteady  one -dimensional  laminar  flame 
propagation  through  a fuel-lean  premixture  under 
a direct  one-step  Irreversible  exothermic  chemi- 
cal reaction,  seem  worth  notinq. 

First,  for  a spatially  unbounded  domain,  the 
effect  of  heat  extraction  from,  or  of  heat  addi- 
tion to,  a flame,  on  the  rate  of  recovery  both  of 
steady  propagation  and  of  flame  temperature  are 
recalled  to  be  as  follows.  The  rate  of  recovery 
is  faster  for  larger  Lewis -Semenov  number,  and 
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Fig.  10c -From  integral -method  results  for  an 
adiabatic  wall,  profiles  of  the  dimensionless 
fuel  consumption  at  the  thin  flame,  -3Y(y,t)/9t 
as  a function  of  tine  t,  for  several  values  of 
the  Lewis-Senenov  number  Le.  Here  3 = 10,  \v  = 1, 
K = 6,  yg  = 10.  The  consumption  rate  at  a fixed 
time  increases  monotonically  with  Lewis  number 
after  flame-wall  interaction  begins 


for  flame- temperature  augmentation  (as  opposed 
to  flame- temperature  depression).  For  too  large 
a flame- temperature  depression,  steady  flame 
propagation  is  not  recovered  at  all,  at  least  not 
for  that  time  span  of  practical  interest  for  many 
purposes.  These  results  are  obtained  here  by 
approximate,  integral -method  treatment  of  a thin- 
flame  formulation. 

Second,  also  for  a spatially  unbounded 
domain,  the  transition  between  steady  flame 
propagation  at  one  stoichiometric  ratio  through 
a continuous  change  to  steady  flame  propagation 
at  an  augmented  (or  decremented)  stoichiometric 
ratio,  is  examined  by  means  of  numerical  solution 
of  simultaneous  nonlinear  Volterra  Integral 
equations,  obtained  by  application  of  the  Fourier 
transform  to  the  thin-flame  formulation  (Figs.  2 
through  4).  The  time  required  for  transition 
through  a given  fuel  stratification  decreases 
with  an  Increase  in  either  the  Lewis-Semenov 
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Fig.  10d -From  integral-method  results  for  an 
adiabatic  wall , profiles  of  the  temperature  at 
the  thin  flame  H as  a function  of  time  t,  for 
several  values  of  the  Lewis-Semenov  number  Le. 
Here,  6 =*  10,  vF  = 1 , K = 6,  Y0  = 10.  The  flame 
temperature  at  a fixed  time  increases  monotoni- 
cally  with  Lewis  number  after  flame-wall  inter- 
action begins.  Since  R(t)  = 1 for  Le  = 1,  the 
parameters  vF,6play  no  role  in  this  case 


Fig.  11  - From  integral -method  results  for  an 
adiabatic  wall,  profiles  of  the  (normalized, 
stoichiometrical ly  adjusted)  fuel  mass  fraction 
Y as  a function  of  the  coordinate  o(=  ip/y ) , at 
time  t = 7.52,  where  Yn  = Y(t=0)  = 10,  for 
several  values  of  the  Lewis-Semenov  number  Le. 
There  is  a slight  depression  of  the  value  of  Y 
for  Le  = 1;  i.e.,  there  is  a nonmonotonic  varia- 
tion (according  to  the  approximate  method) 


number,  or  the  reaction  order  with  respect  to 
fuel,  or  the  dimensionless  Arrhenius  activation 
energy,  or  the  sharpness  of  the  fuel  stratifica- 
tion. Also,  the  time  requires  for  transition  is 
less  for  propagation  through  a fuel  enrichment  of 
a given  magnitude,  as  opposed  to  propagation 
through  a fuel  decrement  of  the  same  magnitude. 
Further,  judicious  use  of  a steady-flame-propa- 
gation relation  yields  the  flame  speed  through 
the  stratification,  with  error  on  the  order  of  a 
few  percent;  that  is,  a quasisteady  approximation 
normally  should  suffice  to  describe  the  flame 
speed  and  flame  temperature  through  the  transi- 
tion. Incidentally,  with  knowledge  of  the  flame 
speed  and  flame  temperature  as  a function  of 
time,  should  the  entire  spatial  profile  for  fuel 
and/or  temperature  be  required,  it  seems  numeri- 
cally easier  to  return  to  the  differential- 
equation  formulation  to  obtain  results,  rather 
than  to  seek  such  information  from  the  integral- 
equation  formulation. 

Third,  attention  is  turned  to  results  from 
integral -method  treatment  of  the  retardation  of 
flame  propagation  as  a cold  wall  is  approached 
(Figs.  5-8).  The  cold-wall  problem  (and  other 
problems)  have  already  been  treated  by  Buck- 
master  [31]  by  the  thin-flame  approximation. 


However,  we  are  unable  to  infer  from  his  paper 
the  results  in  any  of  our  figures  5a-5e,  6a-6c, 
7,  and  8,  except  that  for  the  flame  position  y 
at  Lewis-Semenov  number  Le  = 1 in  Fig.  6a.  We 
have  included  the  details  of  our  treatment  be- 
cause we  believe  the  differences  in  methodology 
to  be  of  interest.  The  results  in  our  figures 
are  more  readily  compared  with  solutions  fur- 
nished by  Kurkov  and  Mirsky  [10]  and  by  Adamczyk 
and  Lavoie  [14],  who  use  numerical  integration 
of  the  initial/boundary-value  problem  with  a 
conventional  law-of-mass-action  expression  for 
the  chemical -reaction  term.  The  position  of  the 
propagating  flame  Y appears  to  reach  a position 
of  closest  approach  to  the  wall  and  then  to 
retreat;  the  reversal  occurs  further  from  the 
wall  at  larger  values  of  the  dimensionless 
Arrhenius  activation  energy  e.  Because  the 
larger-acti vation-energy  flame  maintains  a 
greater  distance  from  the  wall,  the  flame 
temperature  H decreases  less  rapidly  in  time. 
Raising  the  reaction  order  with  respect  to  fuel 
Vr  also  leads  to  closer  approach  of  the  flame  to 
tne  wall  and,  hence,  to  a more  rapid  decrease  in 
flame  temperature.  Increasing  the  Lewis-Semenov 
number  Le  likewise  leads  to  more  rapid  decrease 
of  the  flame  temperature  owing  to  greater  heat 
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transfer  to  the  cold  wall.  Two  other  properties 
of  the  results  are:  (1)  the  rate  of  reactant 
consumption  at  the  thin  flame,  -9Y('f,t)/»t, 
decreases  inonotonical ly  in  time  as  the  cold  wall 
is  approached;  and  (2)  the  speed  of  the  flame  in 
laboratory  coordinates  increases  significantly 
with  flame  exothermicity  K in  the  bulk  gas,  but 
as  the  wall  is  approached  and  chemical  activity 
is  reduced,  the  distinction  in  the  flame  speed 
owing  variation  in  K is  much  less  prominent. 

Fourth,  from  i ntegral -method  treatment  of 
flame  propagation  toward  an  adiabatic  wall  (Figs. 
9-11),  there  is  a terminal  acceleration  of  the 
flame  toward  the  wall,  increasingly  pronounced 
as  the  Lewis-Semenov  number  Le  achieves  values 
in  excess  of  unity;  the  total  amount  of  fuel  is 
consumed  in  finite  time.  The  preheating  of 
remaining  fuel  is  greater  for  larger  Lewis- 
Semenov  number,  and  the  flame  temperature  can 
exceed  the  adiabatic  flame  temperature  for 
Le  > 1;  for  Le  < 1,  the  flame  temperature  falls 
below  the  adiabatic  flame  temperature  as  the 
wall  is  approached.  The  wall  temperature 
approaches  the  thin-flame  temperature  as  the 
flame  reaches  the  wall. 

There  are  two  questions  concerning  accuracy 
that  arise  in  connection  with  the  integral -method 
treatment  of  the  thin-flame  formulation  of  flame/ 
wall  interaction:  first,  how  wel 1 does  the 
approximate,  integral -method  solution  satisfy  the 
thin-flame  iniitial/boundary-value  problem;  next, 
how  wel 1 does  the  approximate,  integral -method 
solution  satisfy  the  initial/boundary-value 
problem  with  the  law  of  mass  action  and  the 
Arrhenius  expression  for  the  specific  rate 
constant.  The  error  is  defined  here  as  the  ab- 
solute value  of  the  residual  obtained  by  substi- 
tuting the  integral -method  solution  into  the 
initial/boundary-value  problem,  normalized  by 
the  absolute  magnitude  of  the  largest  term.  For 
the  Landau  (flame-stabilized)  coordinate  o,  there 
are  three  terms  (temporal  variation,  convective 
transport,  diffusive  transport)  for  the  thin- 
flame  formulation;  there  is  a fourth  term 
(chemical  reaction)  for  the  law-of-mass-action 
formulation.  An  error  may  be  defined  at  any 
spatial  point  in  the  flow  domain  0 < o < ® at 
any  time  t > 0.  (The  accuracy  of  the  integral - 
method  flame/wall-interaction  results  may  be 
checked  also  be  means  of  the  Vol terra- integral - 
equation  formulation,  presented  above  both  for  a 
cold,  isothermal  wall  and  also  for  an  adiabatic 
wall;  this  procedure  for  checking  accuracy  is 
not  adopted  here.) 

The  error  characterizing  how  well  the 
integral -method  solution  approximates  the  solu- 
tion to  the  thin-flame  formulation  seems 
acceptably  small.  For  the  preponderance  of  the 
flow  domain  in  space  and  time,  the  error  is  5% 
or  less.  At  long  times  for  a = 0.9,  the  error 
increases  to  about  20%,  but  this  is  a very 
localized,  peak  value.  For  most  of  the  spatial 
and  temporal  domain,  the  error  characterizing 


how  well  the  i ntegral -method  solution  approximates 
the  solution  to  the  formulation  with  the  law  of 
mass  action  is  approximately  the  same  as  the 
error  incurred  for  the  thin-flame  formulation; 

i.e.,  the  error  remains  small  because  the  reac- 
tion-rate term  is  small.  However,  for  o = 0.8, 

0.9  at  large  times,  the  error  can  become  quite 
large  [i.e.,  0(40:)];  this  is  especially  true  for 
a cold  wall,  since  the  quenched  flame  tends  to 
be  spatially  diffuse.  Nevertheless,  the  error 
incurred  is  quite  localized. 

Worthwhile  directions  for  related  further 
theoretical  and  experimental  work  include  iso- 
baric,  spatially  two-dimensional  flame/wa1’ 
interaction,  and  nonisobaric  effects  on  f.  e 
propagation  through  premixtures  confined  by 
variable- volume  containers.  The  flame-wi thout- 
structure  model  seems  a tractable  means  of  exam- 
ining such  aerothermachemical  flows.  The 
application  of  transform  techniques  for  the 
solution  of  such  thin-flame  models  appears  to  be 
an  approach  deserving  consideration. 
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